Quantum Data Bus in Dipolar Coupled Nuclear Spin Qubits 

Jingfu Zhang/ Michael Ditty/ Daniel Burgarth/ Colm A. Ryan/ C. M. Chandrashekar/'^ 
Martin Laforest/ Osama Moussa/ Jonathan Baugh/ and Raymond Laflamme^'^ 
^Institute for Quantum Computing and Department of Physics, 
University of Waterloo, Waterloo, Ontario, Canada N2L 3C1 
^IMS and QOLS, Imperial College, London SW7 2BK, UK 
'^Perimeter Institute for Theoretical Physics, Waterloo, Ontario, N2J 2W9, Canada 

(Dated: July 28, 2009) 

Abstract 

We implement an iterative quantum state transfer exploiting the natural dipolar couplings in a 
spin chain of a liquid crystal NMR system. During each iteration a finite part of the amplitude 
of the state is transferred and by applying an external operation on only the last two spins the 
transferred state is made to accumulate on the spin at the end point. The transfer fidelity reaches 
one asymptotically through increasing the number of iterations. We also implement the inverted 
version of the scheme which can transfer an arbitrary state from the end point to any other position 
of the chain and entangle any pair of spins in the chain, acting as a full quantum data bus. 

PACS numbers: 03.67.Lx 
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I. INTRODUCTION 



In quantum computation and quantum communication the transfer of an arbitrary quan- 
tum state from one qubit to another is a fundamental element. The most obvious method 
to implement quantum state transfer (QST) on an array of qubits is based on a sequence of 
SWAP gates for neighboring spins. In spin qubit systems the SWAP gate (up to a known 
phase factor) can be implemented through the evolution of the dipolar coupling between 
the neighboring spins for 1/2D time by decoupling the other spins, where D denotes the 
dipolar coupling strength. In experiments, however, the required decoupling operations are 
hard to implement if the spins cannot be individually addressed by spectral selectivity, e.g., 
in large-size solid-state NMR systems. This makes the direct implementation of such gates 
in a large spin system challenging. 

To overcome this problem, schemes based on "always on" spin systems were proposed 
. The state can be transferred with unit fidelity in engineered spin chains or networks 
with XY interactions 3|. However, the required fine-tuned XY couplings are not found in 
natural spin systems In other schemes based on spin chains with Heisenberg interactions 
or with a double-quantum Hamiltonian the fidelity of the QST cannot approach 
unity in scalable systems. 

The above limitations can be relaxed significantly by applying gate operations to receive 
and store the transferred state 6|, |7|. The gates are only applied to two spins at one end of 
a spin chain. In this paper, we experimentally implement the QST in a liquid crystal NMR 
system based on this scheme. Opposed to previous experimental implementations 3] where 
the required XY interactions were engineered by radio- frequency pulses and scalar couplings, 
the dipolar couplings exist naturally in the system and are directly exploited for the QST. 
The dipolar couplings are much stronger (up to 2-3 orders of magnitude) than the scalar 
couplings and therefore can significantly speed up the implementation of the logical gates 
for quantum information processing jsj. The transfer with high fidelity is achieved in an 
iterative manner. Each iteration transfers a finite part of the input amplitude to the target 
spin at the end of the chain. The fidelity of the transfer asymptotically approaches unity by 
increasing the number of iterations. We also experimentally demonstrate the time-inverted 
version of j^. Through this, a full quantum data bus is implemented, where arbitrary 
unknown quantum states can be steered to an any position of the chain. This is also useful 
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for the selective excitation of one spin, which is addressed by the two-spin gates, rather than 
w its individual properties, e.g. chemical shift in NMR. As opposed to previous schemes 
9|, global control is not required. Surprisingly the reversal operation can also be used to 

entangle any pair of spins in the chain by operations at its end only. We demonstrate the 

entangling operation in the qubits at the end points of the chain. 

The QST and its reversal operations mean that the chain is really used as a wire with an 

input, an output, and no gates in the middle; the many-body Hamiltonian of the chain is 

responsible for the transport. Only two spins at the end are required to address. The fidelity 

of transfer converges exponentially fast to unity with respect to the number of iterations. 

The required number of iterations to achieve a good fidelity (e.g. larger than 0.999) scales 



roughly linearly with the system size 



^. Moreover this method is stable when the engineered 



Hamiltonian in implementation deviates the required Hamiltonian 1^. Hence our method 
scales favorably with the size of the spin chains and suitable for large-size systems, such 
as solid or liquid crystal NMR systems, where the differences of the chemical shifts are too 
small to address all the spins individually. 



II. ITERATIVE TRANSFER ALGORITHM IN A SPIN CHAIN 

Our first goal is to transfer the state a|0) + (3\1) from spins j to in a A^- spin chain. 
The Hamiltonian for spins 1 to A^ — 1 is represented as 

1 N-l 

^=2^ E D,,{2aia':-aia'^-ala';), (1) 

j,k=l;k>j 



where a^, cx^ and cr^ denote the Pauli matrices with? indicating the affected spin. Noting 
that H preserves the total number of excited spins [l[ Q], we have 



Ur\0) = e'^\0) (2) 



Ur\i) = afclk) (3) 



k=l 



where Ur = e"*"^^ and e*^ is the (1,1) matrix element of Ur- The state |0) denotes all spins 
pointing up, and |j) denotes all spins up except the spin j pointing down. 

The main operation is the two-spin gate applied only on spins A^ — 1 and A^, and in 
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FIG. 1: (Color online) The numerical simulation (solid) and experimental results (data marked 
by "x") for the probability p„ of the QST as a function of iterations in the four spin system 
used in experiments (see text for dipolar couplings) for transferring a state from spins 1 to 4 (a) 
and entangling the two spins (b) when r = 2.1 ms. In figure (b) j>„ can be approximated as the 
observable coherence Cn [dot-dashed, see Eq. (I19p ] where \Cn — Pn\ ^ 0.0175 when n > 2. The 
experimental data can be fitted as 0.65p„ and 0.77C„, shown as the dashed curves, respectively. 



iteration n the gate is denoted as 



W{Cn, dn) = Il,2,.. 



.N-2 ' 



10 

rf; c; 

-Cn dn 

1 



(4) 



N-1,N 

where /i,...,Ar-2 denotes the unit operator for spins 1 to — 2. The basis order for spins 
AT - 1 and AT is |00), |01), |10) and |11). Noting jc^p + |d„|2 = 1, one finds 



iy(c„,rf„)(c„|N-l) + rf„|N)) = |N). 



(5) 



The N spin system is initialized into the input state a|0) + by setting spin j in the 
system to state a|0) Here j is the location of the sender (receiver) of the QST for the 

(inverse) protocol, which is on some arbitrary spin of the quantum data bus. It is sufficient 
to only discuss the transfer of |j) because Ur only introduces a known phase factor before 
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|0) [see Eq. ([2])] and W{cn,dn) does not change |0) [see Eq. @]. Iteration n is represented 

as 

Qj,n = [h,...,N-2<S)WiCn,dn)][Ur<S>^N]. (6) 

After n iterations one obtains 

N 

IV^n)=Tj,„|j) = ^A,„|k) (7) 
fc=l 

using Eqs. I^Mj- Here Tj^„ = Qj,n---Qj,2Qj,i, AN-i,n = 0, and 

AN,n = a/P^ (8) 

where 

p„=p„_i + |(N-l|[/,(g)/;v|^n-i)|' (9) 
with po = and \ipo) = |j). W{cn, dn) is obtained by setting 

dn = e'^v^/v^ (10) 

Cn={'^-l\Ur^lN\i^n-l)IVp~n- (11) 

In strict nearest-neighbour chains it can be shown that p„ converges to unity by increasing 
the number of iterations. In the present case we have also non- nearest neighbour interactions, 
but numerical results show p„ still approaches unity, with a convergence speed which depends 
on the evolution time r. [See Figure [1] (a)]. The process of QST after a large number of 
iterations can be presented as 

Tj,„(«|0) + m) ^ «e^"'|0) + /5|N), (12) 

i.e., spin ends with the state ae*"^|0) + and e*"^ is known. 

We can exploit the inversion of Tj,„ to implement the QST from spin to spin j, i.e., 
without applying the external operation directly on the spin j to evolve it into state a|0) + 
Hence the spin chain functions as a quantum data bus, which can transfer arbitrary 
unknown states to any qubit. This method also allows to create a selective excitation that 
does not require spectral selectivity, e.g. chemical shift in NMR, to address spin j. The 
external operations are only applied to spins A^ — 1 and A^. By taking the inner product of 
Eq. ([7]) with (N| and using Eq. ([8]) one obtains 

p. = |(j|T^l|N)|^ (13) 
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i.e., Pn is the fidelity for generating |j) by applying T^J to |N). The creation of the selective 
excitation for spin j is represented as 

Trj(a|0) + /3|N)) ^ ae-*"^|0) + (14) 



By modifying the input state one can obtain Tjy„^(ae*"^|0) + /9|N)) a\0) + 

The method of inverse QST furthermore can be used to entangle arbitrary spins j, k 
indirectly by acting at spins — 1 and only. This can be done by designing a pulse 
analogously to Eq. ([6]), and the required pulse sequence is very similar to the inverse QST. 
For this purpose, we set the input state as an entangled state of a pair of spins j and k, 
represented as 

= + \k))/V2. (15) 

Iteration n can still be represented as Eq. Qj, where Qj.n is rewritten as Qj,k,n, noting 
that it depends on the input state. W{cn, dn) is obtained in a similar way by changing 
Po = KNlV^jfc)!^ and \ipo) = \ipjk)- After a large number of iterations we obtain 

Tj,k,n|V^,fc) - |N). (16) 

where Tj k „ = Qj,k,n---Qj,k,2Qj,k,i- From Eq. f|T6|) one can entangle spins j and k with high 
fidelity by applying T^k^,n 1^)) represented as 

^j:k,n|N) - 1^,.). (17) 

The fidelity for generating \^jk) is also represented by Eq. f|T3l) through replacing T^^ by 
T^k,ny |j) by I'ipjk)- The numerical simulation for pn is illustrated as Figured] (b). 

III. EXPERIMENTAL RESULTS 

We use the four protons in ortho- chlorobromobenzene (C6H4ClBr) dissolved in the liquid 
crystal solvent ZLI-1132 as four qubits to implement the experiments. The Hamiltonian is 
represented as 

4 -j^ 4 

HNMR = -7lJ2^^< + l^^ E D,ki2aia':-aia',-ayy). (18) 

i=l k=2,j<k 



Through fitting the spectra [12| obtained by Cory48 [13|] and 1-D MREV-8 pulse sequences, 
and referring to the spectra of molecules with similar structures [s 



we measure ui 
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FIG. 2: (Color online) NMR spectra (blue thick) obtained by Cory48 pulse sequence from the 
thermal equilibrium state (a), and by a collective 7r/2 pulse from pini (b). The red thin spectra 
show the results by simulation. The plot's vertical axes have arbitrary units. The NMR peaks 
marked by "+" indicate the single-quantum transitions between magnetic quantum numbers 2 and 
1. 

106.2, z/2 = —187.7, z/3 = —58.6, 1^4 = 91.3 with respect to the transmitter frequency, 
D12 = -1233.7, Di3 = -149.4, = -93.2, D23 = -716.0, D24 = -236.6, Du = -1677.5 



Hz, and the effective transverse relaxation times (Tg*) as 91, 87, 88 and 82 ms 15|. The 



NMR spectrum obtained by Cory48 from the thermal equilibrium state pth = Yli=i is 



1)(1111| 



16| Hd 



shown in Figure [2] (a). 

All experiments start with the deviation density matrix pj„j = |0000) (0000 1 — 
which can be prepared by the double-quantum coherence Hamiltonian 
^'^ J2k=2,j<k '^jki'^x'^x " '^y'^y) ^ molecuIc with symmetry [l^]. However, we choose 
to generate the effective Hd using a GRadient Ascent Pulse Engineering pulse jisl. Using 
temporal averaging, we prepare pini by summing^ the three states UdPthU\, U\pthUd, '^Pth, 
where Ud = e"**'*''^'' by choosing td = 8.00/1^^2 [iTj]. In the numerical simulation we prepare 
Pini with fidelity 99.97%. Figure [2] (b) shows the NMR spectrum obtained by a collective 
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FIG. 3: (Color online) NMR spectra (a-d) for implementing the QST from spins 1 to 4 after 
100 iterations, when the input states are chosen as cr^ 1 000) (000 1, ay|000)(000|, a2|000)(000| and 
/|000)(000| respectively, where the readout operation e^™y^'^ is applied to obtain observable signals 
in (c) or (d). The plot's vertical axes have the same scale. 

7r/2 pulse in the experiment when the system lies in pini- The NMR peaks marked by "+" 
indicate the single- quantum transitions between magnetic quantum numbers 2 and 1. 

We demonstrate the QST by transferring po from spins 1 to 4 by choosing po = o'x, o-y, 
ijz and /, respectively. Because Tj^„ is spin-preserving, the transitions marked by "+" in 
Figure [2] (b) can represent the QST starting with the input state po|000) (000|. We therefore 
can ignore |1111)(1111| in pj„j and omit the negative frequency spectral region. 

The input state is prepared by applying an operation Uini to pini- With increasing n, 
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FIG. 4: (Color online) NMR spectra for implementing the selective excitation and quantum data 
bus for spin 2. 

Ti,„ transforms po|000)(000| to |000)(000|p asymptotically, where p = (>inea.i2 _ 
experiments we removed the phase factor between p and po by phase correction. For a fixed 
n, we implement the unitary Tx^JJini using one GRAPE pulse. The experimental results 
of the QST after 100 iterations for the various input states are shown as Figures [3] (a-d), 
respectively. 

Exploiting the transformation between the computational basis and energy eigenbasis, 
and ignoring the difference of of the four protons, we can approximatively obtain A^^^ in 
Eq. ([7]) through measuring the amplitudes of the peaks marked by "+" in Figures [3] (a-c) by 
choosing the signals in Figure [2] (b) as the reference. Therefore we obtain p„ = |yl4 „p. For 
the input states a^|000) (000|, (Ty|000) (000|, a,, 1 000) (000 1, pioo is measured as 0.654 ± 0.046, 
0.660 ± 0.052 and 0.693 ± 0.037, respectively. All other |Afc,iooP are below 0.02. 

To observe p„ increasing with n, we also implement the QST by choosing various n when 
the input state is cr2:|000)(000|. The measured p„ is shown in Figured] (a) as the data marked 
by " x", which can be fitted as 0.65p„. 

Next we implement the selective excitation / quantum data bus for spin 2. The reverse 
QST starts with the input state |000) (OOOlcr^; obtained by applying i?^ = e~'"'y'^l'^ to pini- 
When n = 100, T^^l transforms |000)(000|a^ to |0)(0|p|00)(00| with probability close to 1, 
where p = e-«'"^''^z/2^^g«n,e(T2/2_ r^j^^ experimental results are shown in Figure HI The fidelity 
of excitation is measured as 0.744 ± 0.036. 

We choose l^/^u) = (|1) + |4))/-\/2 as the target to demonstrate the entangling operation 
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FIG. 5: (Color online) NMR spectra for measuring the fidelity of the generation of (|1) + |4))/v2. 



in spins 1 and 4. To measure the fidelity, we rewrite Eq. ([T3|l as Pn = | (0000|\I/.„) p [19|] by 
replacing |j) by |^/'i4). Here | \l/n) = -P^T'j^l.nl^) where P denotes the operation to prepare 
I^Pm) from |0000) (e.g. see |20|). When pn is close to 1, we can obtain pn approximately 
by applying a readout operation e^'^'^y^'^ to Noting that in pi^i does not 

contribute observable signals for measuring p„, we approximate p„ as the coherence 



Cn = |Tr(2|0000)(1000|p„) 



(19) 



where p„ = Utot,nPiniUtot,n with Utot.n = e*'^°"i'/'^P"''T]^4 „e*'^°"^/^. The simulated and measured 
Cn is shown in Figure [U (b). The experimental data can be fitted as 0.77C„. Figure [5] 
illustrates the NMR spectra when n = 8. 

The operations Ud, U\, Ry, Tx^nUmi, ^z^n-^y, and Utot,n are experimentally implemented 
using the GRAPE pulses with fidelities in theory larger than 0.99, respectively. The pulse 
lengths are 10 ms for Ud and f/], 20 ms for the other pulses. The experimental errors could 
mainly result from the inhomogeneities of the magnetic field, imperfect implementation 
of GRAPE pulses and decoherence. In order to estimate the quality of the experimental 
spectra, we also list the ideal ones in simulation, shown as the red thin curves in Figures 
EE 



IV. CONCLUSION 



We have given an NMR implementation for various important tasks of quantum control 
that in principle can be achieved indirectly by controlling the end of a spin chain. The 
dipolar couplings naturally existing in the liquid crystal NMR system are directly exploited 
for the QST. The experimental results demonstrate the successful control of the spin system 
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with dipolar couplings by the GRAPE pulses. Firstly, we implemented the transfer of an 
arbitrary quantum state. Secondly, by implementing the reverse QST, we have created a 
full quantum data bus which is controlled by the two-qubit end gates. Finally as another 
application of the reverse QST, we proposed and demonstrated a new method to implement 
an entangling operation. 
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